Thermodynamical Properties of a Rotating Ideal Bose Gas 
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In a recent experiment, a Bose- Einstein condensate was trapped in an anharmonic potential which 
is well approximated by a harmonic and a quartic part. The condensate was set into such a fast 
rotation that the centrifugal force in the corotating frame overcompensates the harmonic part in the 
plane perpendicular to the rotation axis. Thus, the resulting trap potential became Mexican-hat 
shaped. We present an analysis for an ideal Bose gas which is confined in such an anharmonic 
rotating trap within a semiclassical approximation where we calculate the critical temperature, 
the condensate fraction, and the heat capacity. In particular, we examine in detail how these 
thermodynamical quantities depend on the rotation frequency. 
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I. INTRODUCTION 



The rotation of a quantum fluid leads to many inter- 
esting problems. In particular, the rotation of ultracold 
Bose gases has motivated a lot of scientific work within 
the last few years. 

The rotation has the effect of adding an angular momen- 
tum to the system so that it can be compared with the 
motion of a charged particle in a magnetic field. Initial 
experiments [1, 2] have shown that, analogously to the 
superfluid helium II [3] , vortices nucleate at some critical 
rotation frequency. The main difference between earlier 
studies about quantum fluids and EEC's is the trap in 
which the condensate is confined. This inspired A.L. Fet- 
ter [4] to suggest adding a quartic term to the harmonic 
trap potential. His idea was to rotate the condensate so 
fast that the centrifugal force may overcompensate the 
harmonic trapping potential. In a harmonic trap, the 
condensate gets lost when the rotation frequency comes 
close to the harmonic trap frequency [5], but the addi- 
tional anharmonicity ensures the confining of the conden- 
sate. Such a trap was realized in an experiment in Paris 
at the Ecole Normale Supericure (ENS) in the group of 
J. Dalibard in 2004 [6]. 

The experimental setup was a usual magneto-optical 
trap with the frequencies uj^ = tiJy = 2Tr x 75.5 Hz 
and LUz = 27r X II.O Hz. The additional quartic an- 
harmonicity was generated with a Gaussian laser beam 
propagating in z-direction which created a potential 
^ /w^) with the perpendicular radius 



U = Uacxp{-2r 



' x^ 



y'^ , the laser's waste w — 2^ /im, and the 
intensity Uq ~ {ks x 90) nK. Due to the experimen- 
tally realized condition r±_ < w/2, this potential is well 
approximated hy U = Uq ~ {2Uo/w^)rj_ + {2UQ/w^)r^^. 



In Ref. [6] an amount of 3 x 10^ atoms of ®^Rb is set 
into rotation with another laser beam acting as a stirrer. 
The laser creates an anisotropic potential in the xy-plane 
which rotates with frequency fl. In the corotating frame, 
the resulting trapping potential can be written as 



K.ot(x,r!) = 



M 



-n' 



M 



' ± 



2 2 

uj„z ■ 



(1) 



and 
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where uj± = uj^ - 4:UoM~^^w^-^ = 27r x 64.8 Hz 
M is the atomic mass. The last term in (I) corre- 
sponds to the quartic anharmonicity with k = SUqW^'^ = 
2.6 X 10"" Jm-4. 

In the following we treat the Bose gas in the anharmonic 
trap (1) within the grand-canonical ensemble and deter- 
mine the critical temperature, the condensate fraction, 
and the heat capacity of the Bose gas within a semiclas- 
sical approximation. In our discussion the rotation fre- 
quency Q appears as a control parameter. The Paris ex- 
periment [6] allows rotation frequencies f2 up to I.04x w^. 
However, in the present theoretical discussion, we con- 
sider arbitrarily large rotation frequencies. In Fig. I we 
depict how the trapping potential (1) varies with increas- 
ing rotation frequency Q. For small rotation frequencies 
fl < Lu_i_, the potential (I) is convex, for the critical rota- 
tion frequency fl = u!± it is purely quartic in the perpen- 
dicular plane, and for a fast rotation frequency Q > uj± 
the trap has the shape of a Mexican hat. In our dis- 
cussion, we focus on two particular rotation frequencies, 
namely the critical rotation frequency f2 = w^ and the 
limit of an infinite fast rotation frequency fi — > 00 as 
the potential reduces to power laws in both cases. Such 
potentials were investigated some time ago as they lead 
to analytic formulas for the respective thermodynamical 
properties [7, 8]. Note that the case fi ^ cxi corresponds 
to a trap where the bosons are confine d to a cylinder of 
radius Tcyi = ^yM'^uJz{^'^ — uj'j_)/{kh). Thus, the above 
mentioned experimental restriction rcyi < w/2 allows to 
determine a maximum rotation frequency l^max for which 
the anharmonic potential (1) of the Paris experiment is 
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FIG. 1; Trapping potential (1) in (x, 0, 0)-direction 
for the values of the Paris experiment [6] and 
varying rotation frequencies fii = 0, f22 = (^±, 
^3 = 1.04 X (jJx- The latter is the largest exper- 
imentally realized rotation frequency. 



valid. The resulting value f^max ~ 1-08 x lo\_ shows that 
considering an infinite rotation frequency f2 ^ 00 is not 
suitable for this experiment. 



II. IDEAL BOSE GAS IN ROTATING TRAP 

We consider TV particles of an ideal Bose gas vifhich are 
distributed over various quantum states v of the system. 
These states are characterized through the population 
n-n of the one-particle state n of the trap (1), such that 
the energy levels are Ei^ = X^n""-^"' where En denotes 
the one-particle energy. Correspondingly, the number 
of particles in state v are given by N^, ^ X^n^"- "^^^ 
resulting grand-canonical ensemble is determined by its 
partition function 



Z = ^cxph/3(i?,-MiV,)], 



(2) 



where (3 — l/lksT) denotes an inverse temperature, 
ks is Boltzmann's constant, and /i is the chemical po- 
tential. The corresponding grand-canonical free energy 
!F = —{1/(3) \nZ allows to calculate all relevant thermo- 
dynamical quantities [9, 10]. Around the minimum, the 
trap (1) has a small curvature so that the energy lev- 
els are close to each other. With increasing rotation fre- 
quency, the curvature decreases until the critical rotation 
frequency f2 = io± is reached. Not until the rotation fre- 
quency overcompensates the harmonic part of the trap, 
the curvature increases again. Thus, for all experimen- 
tally realized rotation frequencies < J7 < 1.04 x uj±, 
our system can be described by the discrete ground state 
Eq, which must be retained quantum-mechanically, plus 
a continuum of states above Eq. Within this semiclas- 
sical approximation, we can set the ground-state energy 
Eq to zero so that the grand-canonical free energy of the 
ideal Bose gas reads 




FIG. 2: Critical temperature versus rotation fre- 
quency Q. The highest lying line (solid) corre- 
sponds to the data of the Paris experiment [6] , see 
data below Eq. (1). For the deeper lying lines we 
varied the anharmonicity parameter k: k ^> fe/10 
(long dashes), k — * fc/100 (short dashes), and the 
harmonic limit fc | (gray solid). 



^^exp{-/3j[iJ(x,p)-Ai]}, (3) 



00 ^ 

~j-{TjJ (27r?i) 

where the energy levels are replaced by the classical 
Hamiltonian 



^(x,p) = ^ + T4ot(x,f!). 



(4) 



The critical chemical potential ^c, where the condensa- 
tion emerges, is determined by the condition i?(x, p) — 
/ic > and is therefore given by ^c — niiux VJ-ot (x, 51) . 
Due to (1) it reads explicitly 










(-i 



n' 



; rj < uj±_ . 



Q > LUj_ . 



(5) 



Performing the phase-space integral in (3), we obtain 

.F-iVo(Mc-A^) l^i^HHl^ , (6) 

where we have introduced the generalized ^-function 






n' 



X exp 



JlT {^ 



n' 



erfc 



Q' 



with the complementary error function 



erfc(z) = —= 



dte 



(7) 



(8) 



J" 



No{fJ.c - IJ) 



We have also shortened the notation with using 7^ — 
M"^ /{AkksT) as another inverse temperature. We re- 
mark that in the limit of a vanishing anharmonicity fc J, 
with an undcrcritical rotation frequency 17 < uji_ the 
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FIG. 3: Condensate fraction versus reduced tem- 
perature. The solid line corresponds to the con- 
densate fraction (19) of a Bose gas in the trap (1) 
for the rotation frequency f2 = and the param- 
eters of the Paris experiment [6]. The dashed line 
corresponds to the condensate fraction at the criti- 
cal Q. — u>± and at the infinite fast f2 ^ oo rotation 
frequency given by Eq. (18). 
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FIG. 4: Condensate fraction versus reduced tem- 
perature. The condensate fraction (19) is evalu- 
ated for various rotation frequencies Q, the solid 
line corresponds to f2 = to which the tempera- 
ture is normalized. The other lines are for Q — 
uj±/\/^ (short dashes), Q — uj± (longer dashes), 
and Q — •v/3/2 x ojx (long dashes). 



generalized ^-function (7) reduces to the polylogarithmic 
function 



oo 



lim C.(z,0) = C.W = 5]-, (9) 



which is related to the Riemann (^-function via 

c.(i) = CM = E ^ . 



(10) 



Furthermore, we note that in the hniit of the critical 
rotation frequency Q -^ uj^, the generahzed C- function 
(7) reads 

o^™ 2 'o2 =\/^n^C-l/2(g) (11) 

and in the limit of an infinite fast rotation it is approxi- 
mated by 






2V^7?E -;j^Tr, ;f^-oo, (12) 



3 



1^-1/2 



III. CONDENSATE DENSITY 

From the grand-canonical free energy (6) we read off 
that the number of particles A^ = —{d!F/dii)T.v of an 
ideal Bose gas is given by a sum of A^ = Nq + Ne of 
particles in the ground state Nq and particles in excited 
states N„: 



N = No + 



C3(e''^^) 



p^n^uj, {ujI - n^) 



(13) 



The critical temperature Tc at which the condensation 
emerges can be found from Eq. (13) by setting Nq = 



and jjL = He- For undercritical rotation frequencies fi < 
uj_\_ and vanishing anharmonicity fc J, 0, we apply (9) so 
that the critical temperature reads 



Tr.= 



flbJ, 



(- 



9?)N 



^K(3) 



1/3 



0. 



(14) 



For a non- vanishing anharmonicity k, the critical temper- 
ature can not be determined explicitly because it appears 
in Eq. (13) transcendentally. However, there are two spe- 
cial cases in which we obtain an analytical expression for 
the critical temperature. At first, for the critical rotation 
il — ijj^, we find with (11) [11]: 



T, 



ks 



AkTi 



1/5 



U(5/2)j 



2/5 



(15) 



Secondly, the limit of an infinite fast rotation frequency 
il ^ oo leads with (12) to the critical temperature 



Tr. - 



Two, 



kh 



ks \nM^LLi' 



2,, ,3 



1/5 



N 



C(5/2) 



2/5 



(16) 



which is by a factor (1/4)^/^ « 0.76 smaller than the pre- 
vious one. A numerical evaluation of the critical temper- 
ature obtained from (13) is shown in Fig. 2 for the values 
of the Paris experiment [12]. For the non-rotating trap, 
we see that the anharmonicity only slightly affects the 
critical temperature. With increasing rotation frequency 
il, the critical temperature decreases and the difference 
between the harmonic and the anharmonic trap is clearly 
seen. At the critical rotation frequency Q = w_l , the crit- 
ical temperature is Tc — 63.5 nK which is about three 
times smaller than the one estimated for the Paris ex- 
periment [6]. 

From the number of particles (13), we also obtain the 
condensate fraction in the temperature regime T < Tc- 
Here the chemical potential coincides with the critical 



one given by Eq. (5) . For undercritical rotation frequen- 
cies n < u)± and vanishing anharmonicity /c J, 0, we use 
(9) so that the condensate fraction is given by 



N 



; fc = 0. 



(17) 



Furthermore, applying (11) and (12) in the cases Q — 
Lu±_ and il ^ oo, respectively, yields with the critical 
chemical potential (5) the following condensate fraction: 



iVo 

N 



= 1- 



5/2 



(18) 



Here, Tc is given by (14) and (15), (16), respectively. In 
general, the condensate fraction is given by 



TVo 

N 



C3(e^^-,r!) 
C3(e'3=^sr!)' 



(19) 



where Pc = ^/{kBTc). In the low-temperature limit 
rj,0, the condensate fraction shows a power-law behav- 
ior which is different in the two regimes of undercritical 
rotation fi < ll!±_ and overcritical rotation il > uj^. Due 
to (9) and (12), we obtain 



No 

N 



fc|C(3) 



u5/2 



T^ ;n<uj±, 



^_MV^fcrC(5/2)^5/2.^ 



(20) 



Ny/kh 



> UJ± . 



The temperature dependence of the condensate fraction 
Nq/N following from (19) is shown in Fig. 3 and Fig. 4. 
From this we read off that the temperature dependence of 
the condensate fraction depends crucially on the rotation 
frequency and is thus not universal. 



IV. HEAT CAPACITY 

The heat capacity follows from the grand-canonical 
free energy !F = U — TS — fiN , where U is the inter- 
nal energy and S is the entropy, according to 



dU 



(21) 



N.V 



Within the grand-canonical ensemble, the heat capacity 
has to be determined separately in the two regimes T > 
Tc and T < Tr. 



A. Gase Phase 

At first, we treat the gas phase where Nq = and 
determine the entropy via the thermodynamical relation 
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FIG. 5; Approach of the heat capacity (24), solid 
lines, and its approximation (28), dashed lines, to 
the Dulong-Petit law, horizontal line. 



S^-{dT/dT)vy. 

S> 7C4(e''^0) 



keN 2C3(e^^,r!) 



.^ ,2 r.2.2 C3(e^^) 



C3(e^^fi) 



(22) 



The internal energy U = T + TS + fiN then follows from 
(6), (13), (22) and reads 



U> ^ 5C4(e^^^) 
NkeT 2(3(6/3'^, r!) 

+7T(coi - n^f 



C3(e^^^!) 



1 



(23) 



Finally, the heat capacity (21) for temperatures above Tc 
is given by 

35C4(e^^^!) 25C3(e'3^f7) 



C> 

ksN 4C3(e^^,f7) 4(2(6/3^, f^) 



2\2 



+7t(c^1-0^) 



llC3(e''^) C2(e^^) 

2Uef^^,n) C2(e'5^,f^) 



2, 2 o2^4C2( e^^)C2(e^M7)-Cl(e^^) 
C2(e'3^^!)C3(e'3^f}) 



+jU^I - ny 



(24) 



To obtain this result, we have determined the derivative 
[d (3 n / dT) N ,v from (13). In the limit fc J, of a harmonic 
trap, the heat capacity (24) reduces to the well-known 
result 



C> ^, C4(e^^^) g C3(e^^l7) 

ksN - ''c3(e'5^^) ^c2(e'3^f7) ' "■ ^ '^ 

Again, both cases O = wj^ and 17 ^ oo yield the same 
analytic expression: 

C> _ 35 ''7/2 [e"'' ^^'^ > ) 
k^ ~ 



1 C5/2 (ePi^+-<^«-^'^'y 
25C5/2(e''-+--(-i-"^)^) 



C3/2 (e 



/3^+7T(t^l-02)2 



)■ 



(26) 
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FIG. 6; Heat capacity versus temperature without 
rotation, reduced to the critical temperature (14) 
of the trap (1) with A; = 0. The dashed line cor- 
responds to the harmonic trap heat capacity (25), 
(34) for the Paris trap (1) with fc = 0. The black 
solid line is the heat capacity (24), (33) for the an- 
harmonic Paris trap [6]. The horizontal lines cor- 
respond to the Dulong-Petit law: harmonic trap 
(dashed) and anharmonic trap (gray). 



Now, we investigate the heat capacity (24) in the high 
temperature limit T — > oo. To this end we use a large 
T-expansion of the generahzed ^-function (7): 

C. {z, n) « e''^ [V^Ff^iojl - n^) - 27T(^i - n^ + ...] 



o2/3m r 



2"^ 



n' 



47T(t^j 



Q 



2n2 



(27) 



Inserting the expansion (27) into the number of particles 
(13), we find for the first order of the fugacity e^^ w 

2NyTKj^/y/TTM^k%T^, so that the heat capacity (24) 
behaves like 

2 4 



C> , 

keN' 



(cul-n'')-- 



32 7r3/2 



-{c^l-n^. (28) 



Thus, the heat capacity approaches the Dulong-Petit law 
in an anharmonic trap liniT-^oo C^ /(ksN) = 5/2 which 
is 1/2 smaller than the corresponding one in the harmonic 
trap. Furthermore, the first fi dependent term in (28) 
changes its behavior, from being larger (O < uj\_) than 
the limit to being smaller {Vt > uj±), see Fig. 5. 
At the critical point, the harmonic heat capacity (25) 
reduces for small rotation frequencies fl < lu± to 



lim 



C> 



TiT, ksN 
In both limits il 



12 



C(4) gC(3) 



C(3) C(2) 

ll!±_ and Q - 
(26) at the critical point is given by 

C> 35C(7/2) 25C(5/2) 



:4.23; k^O. (29) 
oo, the heat capacity 



lim 

TIT, ksN 



4 C(5/2) 4 C(3/2) 



Condensate Phase 



4.14. (30) 



Now we turn to the condensate phase T < Tc where 
the chemical potential is given by (5). For the entropy. 



we obtain 
5< 



knN 



+1t{^] 



n 



2^2^C3(e''^=)-C3(e^^^^) 



(31) 



C3(e^=^S^) 
The internal energy below the critical temperature reads 



c/< 



NksT \tJ l2C3(e^=^^0) 

C3(e^^0-e(^x-^)C3(e''^%^) 



-7T(a 



n 



2n2 



where Q denotes the Heaviside function. Thus, we find 
for the heat capacity below the critical temperature with 
help of Eq. (5) 



kuN 



fry f 35Uef>'^^,n) , 2 



Q 



2\2 



llC3(e 



P^J.c^ 



2C3{e^^^^,n) 



5e{uji_ - n) 



C3(e^^=,f2) 



2\4 



+ -f^{iji-n') 



e{n-ujj_) 



C2(e^^0 



e{ujj_ - n) 



C3(e^^^S^) 
C2(e^^=)-C2(e''^^^)" 



In the limit fc i 0, Eq. (33) simplifies to 



c< ^^^m 



^1 ; fc = 0. 



kBN C(3) \Tc 

Thus, at the critical point it has the value 



lim 

TTTe ksN 



12 



C(4) 
'C(3) 



10.80; fc = 0. 



(33) 



(34) 



(35) 



In both limits Q = ll!±_ and il ^ oo, the heat capacity is 
given by 



C< _ 35 C(7/2) /T Y"^/^ 
fc^"Tc(5/2) [t, 

33) at the 
35 C(7/2) 



(36) 



so that the heat capacity (33) at the critical point reduces 
to 

C< 



T™ keN 4 C(5/2) 



7.35. 



(37) 



In the low temperature limit T I 0, the heat capacity (33) 
tends to zero in accordance with the third law of ther- 
modynamics. We note that the low-temperature limit of 
the heat capacity (38) has the same power-law behavior 
as the corresponding one of the condensate fraction (20): 

354C(4) 



kBN 



rpo 



ANffLo^{Lo\-n^) 
AN^/kh^uj^ 



Vl < U!± 



(38) 



> UJ± 
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FIG. 7; Heat capacity versus temperature, reduced 
to the critical temperature TciO. = 0) of the an- 
harmonic trap (1) for varying rotation frequencies. 
The lines correspond to f2 = (solid), f2 = cjx 
(long dashes), and fi = 2lox_ (short dashes). The 
gray solid line corresponds to the Dulong-Petit law, 
the first term of (28). We note that for fi = 2ci;x 
and T > Tc the heat capacity approaches the 
Dulong-Petit limit from below (see Fig. 5). 



Fig. 6 shows the temperature dependence of the heat 
capacity for the values of the Paris experiment without 
rotation [6] . We see that the effect of the anharmonicity 
is rather small. According to Ehrenfest's classification, 
the discontinuity at the critical temperature characterizes 
the phase transition to be of second order. In Fig. 7 
wc show how the heat capacity depends on the rotation 
frequency fl. Here, the rotation has a huge influence on 
the temperature dependence of the heat capacity. 



perature Tc ~ 64 nK at the critical rotation frequency 
fi = 0J± yields kBTc/TiuUz ~ 120 ^ 1 which justifies 
a posteriori our semiclassical treatment. However, the 
value Tc ~ 64 uK is about three times smaller than the 
one estimated in the Paris experiment [6] . This huge dis- 
crepancy could not be explained with the circumstance 
that our semiclassical analysis of the rotating ideal Bose 
gas neglects three important aspects, namely finite-size 
corrections, interactions between the particles, and the 
effect of vortices in the condensate. All three should 
have the effect of lowering the critical temperature. In 
a harmonic trap, it was found numerically that the 
finiteness of the system slightly lowers the critical tem- 
perature [13] which was also shown analytically [14, 15]. 
Furthermore, an additional weak repulsive two-particle 
contact interaction leads to a negative shift of just a 
few percent [16, 17]. Therefore, we conclude that our 
semiclassical findings ask for a revised experimental 
measurement of the critical temperature for a rotating 
Bose gas. 

For the condensate fraction, we can state that the 
low-temperature behavior crucially depends on the rota- 
tion frequency f2. It shows a non- uniform temperature 
dependence which is in between the two power-laws T^' ^ 
and T^. 

The heat capacity of the rotating ideal Bose gas is 
discontinuous at the critical temperature. It tends in 
agreement with the third law of thermodynamics to 
zero in the low temperature limit T | 0, and approaches 
the Dulong-Petit law (28) in the high temperature limit 
T ^ oo from above {Q < lu±) or from below {Q > lu±). 



CONCLUSIONS 



We have determined a semiclassical approximation for 
the critical temperature Tc at which the condensation of 
a rotating ideal Bose gas occurs in the anharmonic trap 
(1). We have found that condensation is possible even 
in the overcritical rotation regime f2 > uj± which is in 
contrast with the harmonic trap where the condensate 
gets lost, when the rotation frequency fl gets close to 
the trap frequency uj±. Our value for the critical tem- 
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